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Abstract — Due to a variety of applications, there is an ever-increasing demand for urea and subsequently its 
production process remains a popular research topic. In the current climate however, studies for solving 
industrial challenges and the search for a more sustainable process design are required. Previous works 
concerning simulation of industrial scale production have been developed, however almost none of them are 
reproducible nor consider urea quality parameters. The severe process operating conditions and the lack of 
biuret information are the main challenges in modeling and simulating such complex process. This paper 
proposes a systematic approach for simulation and validation of the current urea production process. Industrial 
data from the largest operational urea facility in Latin America are used. Simulation is validated against more 
than SO industrial parameters. Deviation of less than 6% is obtained for mass composition and less than 8% for 
other variables considered. This work is a key point for retrofit studies and design of new processes models. 

Keywords— industrial data validation, industrial process simulation, urea process. 


A mathematical model for synthesis section was 
developed by [6]. The thermodynamic framework was 
based on Wilson and ideal gas equations. Inlet and outlet 
temperatures and mass fractions in the reactor and 
Scrubber outlet were compared to industrial data and 
varied from -6.9% to 2.6%. Zhang et al. (2005) simulated 
also the high synthesis loop. Extended electrolytic 
UNIQUAC equation and perturbed-hard-sphere were 
employed for thermodynamic modeling. Reactor and 
stripper outlet mass fractions varied from 2.7% to 9.7% 
when compared to industrial data. Rasheed (2011) 


There are a number of reported studies ([2]; [3]; [4]; 
[5]; [6]; [7]; [8]; [9]; [10]; [11] and [12]) on mathematical 
modeling and/or simulation synthesis section - the 
reaction section. Nevertheless, there is still a range of 
restrictions for simulating such a complex process. 


Demand for urea is constantly increasing. Widely used as 
nitrogen-based fertilizer, additive in animal feed and in 
cosmetic industries, urea has recently taken a leading role 
reducing NOx emissions for diesel engines [1]. 


I. INTRODUCTION 


simulated the urea reactor applying SR-POLAR equation 
for thermodynamic modelling and proposed a power law 
kinetic for ammonium carbamate and urea formation. 
Deviations from industrial data were reported as less than 
5.0% for liquid composition in the reactor outlet. 
Zendehboudi et al. (2014) proposed a mathematical 
model for urea reactor based in a UNIQUAC approach. 
When compared to industrial data, deviation less than 
2.3% for the liquid outlet stream is obtained. Edrisi et al. 
(2016) simulated the entire urea plant using SR-POLAR 
for thermodynamic modeling. Industrial data deviation 
and biuret reaction weer not reported. Chinda et al. (2017) 
simulated the synthesis loop through SR-POLAR basis 
and proposed a power law kinetic model for ammonium 
carbamate, urea and biuret formation. Deviations from 
industrial data were less than 6%. Jeenchay et al. (2018) 
simulated urea process using NRTL for thermodynamic 
approach and no validation was presented. 


The main difficulty in simulating urea process is still 
the availability of physical-chemical data in the range of 
conditions observed along the entire process. An 
important quality parameter for urea as final product, the 
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biuret content, is lacking in available data at relevant 
process conditions. Just Hamidipour et al. (2005), 
Zendehboudi et al. (2014) and Chinda et al. (2017) had 
considered biuret reaction in synthesis section. Besides 
this, further studies using the developed simulation as 
basis are only presented in [12], as an economic analysis 
of the process. A validated simulation is a reliable way to 
identify industrial bottlenecks in the current urea process 
and a key point for studies aiming in promoting 
innovation and technology breakthroughs for industries. 

This paper proposes a systematic approach for simulating 
and validating urea process. For this, industrial data from 
the largest operational urea facility in Latin America and 
biuret reaction are considered. Employed methodology is 
presented in three steps: Step 1 - Industrial data 
collection; Step 2 - Process Simulation, Step 3 — Process 
validation. 

In order to guarantee reproducibility for other urea 
industrial cases, all steps are performed using commercial 
software and the main simulation parameters are 
presented. 

II. METHODOLOGY 

The methodology presented here is hierarchical and is 
composed of three steps. Each step can be used 
independently given that information from previous step 
is available. 

Step 1.1. Industrial Data Collection 
At this step, all industrial data (flows, compositions, 
temperatures and pressures) are collected. The intention 
of this step is to obtain enough information to model the 
process and validate the simulation. Industries usually 
have their own agenda for collecting data concerning 
stream compositions depending on the analytical 
equipment used and the laboratory procedures and 
schedules. 

It is important to collect data from all available 
composition analyzers and flowmeters in the plant in 
order to validate the mass balance of the simulation. In 
order to validate the energy balance, it is necessary to 
collect data from pressure and temperature indicators. It is 
convenient to have these data from points as close as 
possible to the composition analyzers, such that flow rates 
can be estimated where necessary. 

Step 2. Process Simulation 

The steady state simulation proposed for urea process is 
built in AspenPlus®. Ammonium carbamate, urea and 
biuret reactions are considered, given that urea is formed 
only in liquid phase. Industrial data do not consider the 


ammonium carbamate mass fraction. Therefore, it was 
assumed that 99.0% of CO 2 reacts to form ammonium 
carbamate, according to [8] and [13]. 

Thermodynamic modeling is based on SR-POLAR 
equation, recommended for highly non-ideal systems at 
high temperatures and pressures and for both non-polar 
and high polar components, according to [14], [15] and 
[16]. Kinetic equations are taken from [10]. 

Pure component data and binary interaction parameters of 
NH 3 , CO 2 , H 2 O, urea, ammonium carbamate, N 2 , O 2 and 
H 2 are taken from the AspenPlus database. Biuret pure 
component data is obtained from NIST (National his titute 
of Standards and Technology) and DECHEMA 
(Gesellschaft fur Chemische Technik und 
Biotechnologie) database. In terms of vapor pressure and 
binary interaction parameters, biuret is assigned the same 
parameters are urea. This step results in detailed mass and 
energy balance data and the properties of all streams in 
the flowsheet. 

Step 3. Process Validation 

The main objective of this step is to perform the 
validation of the simulation using the data collected in 
Step 1. Thus, it is necessary to process all the industrial 
data in order to evaluate which data can be used to 
validate the simulation, since industrial data may present 
some fluctuation during operation. All the plant data 
collection, performed in Step 1, should be taken at the 
same time or, at least, on the same day. This is a point to 
be highlighted, given that inter-connected industrial 
plants do not operated at steady-state and it is important 
to ensure stable operating points are used in the data 
treatment. In a urea production complex, it is possible for 
example that the ammonia unit is shut-down before the 
urea plant experiences deviations due to upstream process 
e.g. natural gas/residue asphaltic processing. It is less 
important to understand the nature of the up-streams 
deviations, as long as it is possible to identify deviations 
in the given process data in order to rule them out of 
validation process. For this, an analysis with the variation 
coefficient is performed with the capacity data taken each 
4 hours. Coefficient of variation with values less than 
1.5% are considered to represent data that are not varying 
significantly [10] and, therefore, indicate a steady-state 
condition in the process. These two procedures can 
guarantee that data used for validation correspond to a 
stable and continuous operation. Finally, the validation 
process can be performed calculating the deviation 
between industrial experimental data and data obtained 
from the simulation. Thus, this step can be described as: 
(a) from the processed data from Step 1, select only those 
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ones that were taken on the same day. (b) collect 
production capacity data, (c) calculate arithmetic means, 
sample standard deviation and coefficient of variation 
using the production capacity data, (d) eliminate data with 
variation coefficient greater than 1.5%. (e) evaluate the 
deviation between industrial experimental data and data 
obtained from the simulation. 

III. RESULTS AND DISCUSSION 

The process unit analyzed produces 2000 ton/day of urea 
through Stamicarbon technology and can be divided into 
five blocks: synthesis, evaporation, prilling, desorption & 
hydrolysis and recirculation. A simple block diagram of 
the process can be seen in Fig. 1. 



flowsheet and data point collection can be seen in Fig. 2 
and Fig. 3. Mass composition analyzers are indicated in 
blue, flowmeters in green, pressure indicators in yellow 
and temperature indicators in red. In order to facilitate the 
sequence of processing units, TAG order is based on 
section unit (S-Synthesis; E- Evaporation; P-Prilling; 
D&H-Desorption & Hydrolysis; R-Recirculation) and 
flow streams (numerical sequence). 

Step 2. Process simulation 

Process simulation was performed as described in Step 2. 
From AspenPlus model library: urea reactor was modeled 
as a sequence of CSTRs in series; pool condenserusing 
R-Stoic; stripper, scmbber, rectifying column, absorber, 
desorbers and hydrolyzers as RadFrac columns; main heat 
exchangers, condenser and evaporator T-5 were modeled 
as Heat-X; evaporators T-6 and T-7 were modeled as V- 
dmm list of the main equipment and the correspondent 
AspenPlus model library used for simulation can be seen 
in Table 1. 


Fig.l - Simplified block diagram for industrial urea 
production. 

The main equipment in the synthesis section are: pool 
condenser, reactor, scmbber and stripper. This section is 
responsible for ammonium carbamate and urea 
production. After leaving the synthesis section, the liquid 
product from the stripper is sent to recirculation. This 
section is responsible for removing the ammonium 
carbamate present in the solution through its 
decomposition in NH3 and CO2, besides condensing NH3 
and CO2 into ammonium carbamate and recycle them 
back to the synthesis section. The recirculation section 
consists in: rectifying column and its respective heaters, 
condensers and an atmospheric flash tank. The production 
from the recirculation section is fed to the evaporation 
section in order to concentrate up the urea solution before 
it is prilled. This section operates under vacuum, which 
means a big part of water, ammonia and carbon dioxide 
are removed from the solution. The evaporation section 
consists mainly of three evaporators and its respective 
heaters. The last section is called desorption and 
hydrolysis and consists in an adsorber and two desorber 
units. The main function of this section is to recovery 
NH3, CO2 and urea present in the water that comes from 
the Evaporation section. 

Step 1. Industrial data collection 

Industrial data collection was performed as described at 
Step 1. Thus, mass composition analyzers and flowmeters 
were identified in the industrial flowsheet. Further, 
pressure and temperature indicators closest to mass 
composition analyzers were also identified. Process 


Table 1 - Model library from AspenPlus. 


TAG 

Unit 

Model 

TAG 

Unit 

Model 

S-2 

Stripper 

RadFrac 

R-7 

Rectifying 

RadFrac 

S-4 

Pool 

R-Stoic 

R-8 

Condenser 

Heat-X 

S-5 

Reactor 

RCSTR 

R-12 

Absorber 

RadFrac 

S-6 

Scrubber 

RadFrac 

H& 

D-17 

Desorber 1 

RadFrac 

E-24 

Condenser 

Heat-X 

H& 

D-18 

Hydrolyzer 

RadFrac 

E-25 

Evaporator 

V-druin 

H& 

D-19 

Desorber2 

RadFrac 

E-26 

Absorber 

RadFrac 





Step 3. Process Validation 

Data from a total of 270 operational days were provided 
by an industrial urea plant for validating the simulation. 
As described at Step 3(a), in order to guarantee 
consistency to the analysis, all the experimental points 
should be taken in the same operational condition. Only 
32 operational data points (days) met this criterion, (b) 
Production capacity from this data was taken in intervals 
of 4 hours, (c) Arithmetic mean (AM), sample standard 
deviation (SSD) and variation coefficient (VC) were 
calculated, (d) Production capacity data with coefficient 
of variation with values lower than 1.5% were selected. 
Table 2 presents the data used and obtained at this step 
for points with VC lower than 1.5%. The production 
capacity range varied from 86.45% to 98.21%. (e) 
Validation of the simulation was performed calculating 
the difference between industrial and simulated data and 
dividing it per industrial data. 
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Synthesis 


Recirculation 

Hydrolysis and Desorption 


S-1 C0 2 compressor 


R-7 Rectifying column 

H&D-15 Reflux condnser 


S-2 Stripper 


R-8 Pre-condenser 

H&D-16 Level tank for reflux condenser 


S-3 NH 3 pump 


R-9 Ammonium carbamate pump 

H&D-17 Desorber 1 


S-4 Pool Condenser 


R-10 Ammonium carbamate condenser 

H&D-18 Hydrolyser 


S-5 Reactor 


R-11 Ammonium Carbamate tank 

H&D-19 Dessorber II 


S-6 Scrubber 


R-12 Absorber 

H&D-20 Dessorber heat exchanger 




R-13 Circulation cooler 

H&D-21 Hydrolyser feed heat exchanger 




R-14 Circulation pump 

H&D-22 NH 3 -H 2 0 tank 



Fig.2 - Process flowsheet and data point collection for Synthesis, Recirculation and Hydrolysis & Desorption sections. 
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Evaporation 

E-23 Flash tank 

E-24 Condenser flash tank 

E-25 Pre evaporator 

E-26 Atmospheric absorber 

E-27 Pre evaporation condenser 

E-28 Condenser 1 st stage evaporator 

E-29 Urea tank 

E-30 Urea solution pump 

E-31 Ejector 1st stage 


E-32 1st stage evaporator 
E-33 2nd stage evaporator 
E-34 Concentrated urea pump 
E-35 Booster 

E-36 1st condenser from 2nd stage evaporator 

E-37 1st ejector 2nd stage 

E-38 2nd condenser from 2nd stage evaporator 

E-39 2nd ejector 2nd stage 

E-40 After condenser 

E-41 Prilling tower 


Fig.3 - Process flowsheet and data point collection for Evaporation section. 
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Table 2 - Statistical Analyses with capacity planta data. 


Poin 

t 

cc 

AM n VC 
(%) (%) 

Poin 

t 

cc 

AM n VC 

(%) (%) 

A 

8 ^ 4 0.39 0.45 

G 

95 ‘ 6 0.14 03 

O J 

B 

8 ^ 7 0.39 0.45 

H 

95 ? 8 0.29 

C 

87 -° 1.05 1.21 

0 

I 

96 2 6 0.63 ° 5 6 

D 

873 0.07 0.08 

J 

983 0.11 03 

E 

07 ^ 

0.52 0.60 

K 

983 0.49 

F 

87 4 

0.10 0.11 




Equations used for process validation step, can be seen in 
Table A Appendix A. Fig. 4 shows the comparison of 
mass fraction for each component in the outlet of the 
reactor (R) and the stripper (S). As it can be seen there is 
a great accordance between industrial (HMD) and 
simulation (SIM) results. 
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Fig.4 - Comparison between mass fraction composition in 
the outlet of the reactor and the liquid outlet of stripper. 

A total of 37 different process parameters were evaluated, 
among them stream temperature, steam generation, mass 
fraction for CO 2 , NH 3 , urea, H 2 O and biuret, CO 2 
conversion in the reactor and stripper efficiency. For all 
mass fractions evaluated, the deviation between the value 
predicted by the simulation and the real value obtained 
from industrial data were less than 6%, while for steam 
generation and stream temperatures the deviation was less 
than 8%. A selected list of variables and the comparative 
deviations with other similar works are given in Table 3. 
As it can be se seen for reactor and stripper results were 
very close to reported data from literature. 


Table 3 - Absolute average deviation for evaluated points. 


Equipment 

Parameter 

This 

work 

lite¬ 

rature 

Ref. 

Pool 

LP steam flow 

7.56% 

- 

- 

Condenser 

UreaMF 

0.89% 

- 

- 


CO 2 MF 

5.95% 

8.84% 

[ 6 ] 


NH 3 MF 

4.33% 

9.76% 

[24] 


UreaMF 

3.38% 

2.65% 

[5] 

Reactor 

H 2 OMF 

4.38% 

2.71% 

[24] 

Biuret content 

4.96% 

- 

- 



CO 2 conversion 

2.41% 

0.44% 

[3] 


N/C ratio 

4.83% 

6.90% 

[5] 


H 2 O /urea ratio 

6.18% 

- 

- 

Scrubber 

T of the liquid 
outlet 

4.77% 

- 

- 


CO 2 MF 

5.00% 

5.10% 

[24] 


NH 3 MF 

6 . 12 % 

4.14% 

[24] 


urea MF 

2.93% 

0 . 20 % 

[4] 


H 2 OMF 

4.32% 

4.96% 

[24] 


Biuret content 

4.96% 

- 

- 

Stripper 

Efficiency 

2.53% 

- 

- 


N/C 

2.42% 

- 

- 


H/urea 

6.50% 

- 

- 


Steam flow 

4.27% 

- 

- 


T of liq outlet 

6.16% 

- 

- 


Urea production 

3.56% 

- 

- 

Rectfying 

column 

UreaMF 

1.52% 

- 

- 

H 2 OMF 

3.04% 

- 

- 


Biuret content 

3.18% 

- 

- 

Urea 

tank 

UreaMF 

4.17% 

- 

- 

H 2 OMF 

4.83% 

- 

- 


Biuret content 

3.82% 

- 

- 

Final 

H 2 OMF 

6,91% 

- 

- 

product 

Biuret content 

5,73% 

- 

- 

Reflux 

CO 2 MF 

4.34% 

- 

- 

condensate 

tank 

NH 3 MF 

4.89% 

- 

- 

urea MF 

2.58% 

- 

- 


CO 2 MF 

4.31% 

- 

- 

Ammonia 

NH 3 MF 

4.55% 

- 

- 

water tank 

UreaMF 

3.49% 

- 

- 


H 2 OMF 

0.34% 

- 

- 
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IV. CONCLUSION 

A simulation for all sections of urea production is 
developed and validated against more than 30 industrial 
parameters using a total of 32 processed daily operations 
data. Good consistency between simulation results and 
industrial data is presented, being that a deviation of less 
than 6% is obtained for mass composition and less than 
8% for other variables considered. Reproducibility of 
other industrial urea plants is therefore possible and 
permits using it for reliable retrofit studies and design of 
new processes models 
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APPENDIX A 


Table A - Equations used in validation step. 


Simple Arithmetic 
Average 

x = Y J x l 

Sample Standard 
Deviation 

\n-l , = „ 

Coefficient Of 

Variation 

Cv = — 

X 

Deviation between 


simulated and 

industrial data 

, , _, industrialdata— simulationdata , 

deviationf%) =- ; - ; - x 100 

industrialdata 

CO 2 conversion 

urea 

% c ° 2 u rea + 1,3 65 x C0 2 

N/C relation 

NH 3 2x urea+ 3,53 x NH 3 

CO 2 urea +1,365 x C0 2 

H/U relation 

H,0 H 2 Ox 3,33 

urea urea 

Stripper efficiency 

urea 

^ urea+ 1,765 x NH 3 

Biuret content 

biuretcontent biuretmass [&g] 

biuretmass + ureamass [/eg] 
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